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PREFACE 


Th'1  work  described  in  this  report,  was  performed  under  Contract  No. 

DACW  -\>-i  '-0-00.'  t,  "Cracking  of  Earth  Dams,"  between  the  U.  S.  Army  Engineer 
Waterways  Experiment  Station  (WES)  and  Harvard  University.  The  contract  was 
sponsored  by  the  Office,  Chief  of  Engineers,  under  Engineering  Studies  Item 
ES-SM ,  "Cracking  of  Earth  Pains." 

The  general  objective  of  this  research,  which  began  in  1968,  was  to 
investigate  by  the  finite  clement  method  the  factors  that  influence  cracking 
in  eartli  dams.  Work  on  this  project  \vas  conducted  under  the  supervision  of 
Arthur  Casagrande,  Professor  of  Coil  Mechanics  and  Foundation  Engineering. 
The  project  was  administered  by  the  President  and  Fellows  of  Harvard 
University.  The  report  was  prepared  by  Arthur  Casagrande  and  S.  W. 
Covarrubias . 

The  contract  was  monitored  at  WES  by  Mr.  J.  B.  Palmerton,  Rock 
Mechanics  Section,  Soil  and  Rock  Mechanics  Branch,  Soils  Division.  Mr. 

J.  P.  Sale  was  Chief  of  the  Soils  Division  during  the  preparation  and 
publication  of  this  report.  Contracting  Officers  were,  successively,  COL 
Levi  A.  Brown,  CE,  and  COL  Ernest  D.  Peixotto,  CE,  Directors  of  WES. 
Technical  Director  of  WES  during  this  time  was  Mr.  F.  R.  Brown. 
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FOREWORD 


The  investigations  reported  herein  were  performed  in  partial  fulfill¬ 
ment  of  Contract  No.  DACW  E9-69-C-0029  between  the  U.  S.  Army  Engineer 
Waterways  Experiment  Station  and  Harvard  University,  dated  26  March  1969, 
and  of  the  modification  of  this  Contract  effective  15  August  1969*  The 
general  purpose  and  scope  of  this  Contract  was  stated  as  follows: 

",  .  .to  investigate  factors  influencing  the  development  of 
cracks  in  earth  dams  using  the  finite  element  method  to  de¬ 
termine  stress  and  strain  distribution  in  earth  dams  for  a 
variety  of  typical  boundary  conditions  (in  particular,  various 
shapes  of  abutments),  and  of  stress-strain  properties  of  the 
materials  in  the  dam  and  its  foundations.  Dams  which  have 
cracked  will  be  analyzed  to  establish  empirical  correlations 
between  analytical  results  and  actual  performance." 

This  scope  was  later  enlarged  to  include  the  following  studies : 

(1)  Investigation  of  factors  affecting  the  development  of  tension 
zones  around  conduits  in  earth  embankments . 

(2)  Investigation  of  tension  zones  caused  by  rigid  cutoff  wail  be¬ 
neath  an  earth  dam. 

(3)  Investigation  of  additional  case  records  to  compare  observed 
tensile  strains  along  the  crest  with  the  results  of  analyses 
by  means  of  the  finite  element  method. 

The  original  Contract  was  fulfilled  by  the  doctoral  thesis  of  Sergio 

Covarrubias  entitled  "Cracking  of  Earth  and  Rockfill  Dams;  A  Theoretical 

Investigation  by  Means  of  the  Finite  Element  Method, "  published  as  WES 

Contract  Report  S-69-5,  April  1969-  hems  (l)  and  (2)  above  are  fulfilled 

by  the  present  report  and  item  (3)  will  be  fulfilled  by  a  forthcoming 

report . 
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SYNOPSIS 

The  purpose  of  this  study  was  to  investigate  by  means  of  the  finite  element 
method  (1)  the  effect  of  rigid  conduits  and  cutoff  walls  on  the  stress  distribution  and 
on  the  development  of  tension  zones  in  embankments  and  their  foundations,  and  (2) 
the  distribution  of  stresses  acting  on  the  sides  of  conduits  and  cutoff  walls. 

All  materials  were  assumed  to  be  linearly  elastic  with  equal  properties  n  ten¬ 
sion  and  in  compression.  The  only  load  considered  was  the  weight  of  the  embank¬ 
ment  and  it  was  assumed  to  be  applied  in  a  single  lift. 

Detailed  investigation  by  means  of  the  finite  element  method  included  six 
cases  of  conduits  and  four  cases  of  cutoff  walls.  In  the  case  of  conduits,  it  was 
demonstrated  that  tension  zones  occur  adjacent  to  conduits  with  sharp  edges,  and 
that  the  presence  of  a  zone  of  more  compressible  material  above  the  roof  effectively 
improves  the  loading  conditions  on  the  conduit.  In  the  case  of  cutoff  walls  beneath 
embankments,  it  was  shown  that  tension  zones  develop  in  the  upper  part  of  the  em¬ 
bankment  and  adjacent  to  the  top  of  the  wall.  While  the  presence  of  a  zone  of 
more  compressible  material  on  top  of  the  wall  reduces  the  load  on  the  upper  sur¬ 
face  of  the  wall,  it  somewhat  increases  the  friction  forces  along  the  sides  of  the  wall. 
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CHAPTER  1 
INTRODUCTION 


1.1  PURPOSE  AND  SCOPE 

The  purpose  of  this  study  was  to  investigate  by  means  of  the  finite  clement 
method  (1)  the  effect  of  rigid  conduits  anu  cutoff  walls  on  the  stress  distribu’ion 
and  on  the  development  of  tension  zones  in  embankments  and  their  foundations,  and 
(2)  the  distribution  of  stresses  acting  on  the  sides  of  conduits  and  cutoff  walls. 

Conduits  —  All  conduits  were  assumed  to  be  rigid  and  to  be  supported  direct- 
iy  on  an  incompressible  base  which  also  forms  the  foundation  of  the  embankment. 

The  following  variables  were  investigated:  (1)  shape  of  conduit;  (2)  thickness  of  over¬ 
lying  embankment;  and  (3)  effect  of  a  more  compressible  zone  located  immediately 
over  the  conduit. 

Cutoff  Wails  —  Four  cases  were  analyzed:  (1)  a  thick  cutoff  wall  with  rough 
sides,  i.e.  assuming  no  relative  displacements  between  the  wall  and  the  adjacent  soil; 
(2)  a  thick  slippery  cutoff  wall,  i.e.  assuming  no  shear  stresses  on  the  sides  of  the 
wall;  (3)  a  thick  rough  cutoff  wall  with  a  compressible  zone  on  top  of  the  wall;  and 
(4)  a  thin  cutoff  wall  with  rough  sides.  Particular  attention  was  also  paid  to  the  dis¬ 
tribution  of  normal  and  shear  stresses  along  the  sides  of  the  wail  in  an  effort  to  de¬ 
termine  the  vertical  load  on  the  wall  produced  by  negative  skin  friction. 


U  BASIC  APPROACH 

The  studies  under  Contract  No.  DACW  39-69-C-0029  are  based  on  the 
senior  author’s  concept  that  for  the  purpose  of  investigating  tension  zones  and 
cracking  in  embankments  it  is  of  advantage  to  assume  that  all  materials  are  perfectly 
elastic  and  that  although  such  simplification  of  the  stress-strain  properties  of  the  mate¬ 
rials  exaggerates  the  magnitude  of  the  tensile  stresses,  it  does  not  significantly 
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change  the  geometry  of  the  tension  zones  and  the  locations  of  the  maximum  tensile 
stresses  as  compared  to  those  that  develop  in  an  actual  embankment.  This  hypothesis 
was  thoroughly  investigated  by  the  junior  author  in  his  doctoral  thesis  (Covarrubias, 
1969)*  where  he  compared  observational  data  of  several  dams  with  the  results  of  the 
finite  element  analysis  of  these  dams  and  found  good  agreement,  thus  demonstrating 
the  usefulness  of  this  approach. 


1.3  METHOD  OF  ANALYSIS 

The  method  of  analysis  and  the  computer  program  are  the  same  as  described  by 
Covarrubias  (1969,*  Chapter  3  and  Appendix). 

The  finite  element  method  as  used  in  this  investigation  yields  values  of  stresses 
at  the  centroids  of  the  elements.  Linear  extrapolation  was  used  to  evaiuate  the  stresses 
at  the  boundaries  of  the  sections  analyzed. 


* 


Covarrubias,  S.  W.  (1969),  “Cracking  of  Earth  and  Rockfill  Dams,”  Harvard  Soil  Mechanics 
Series  No.  82,  Harvard  University,  Cambridge,  Mass. 
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CHAPTER  2 

TENSION  ZONES  IN  AN  EMBANKMENT  ON  A  RIGID  FOUNDATION 
CAUSED  BY  A  RIGID  CONDUIT 


2.1  GENERAL 

The  influence  of  the  following  variables  on  the  development  of  tension  zones 
and  on  the  stresses  acting  on  the  conduit  was  investigated:  (1)  shape  of  conduit,  (2) 
thickness  of  overlying  fill,  and  (3)  effect  of  a  more  compressible  zone  located  imme¬ 
diately  over  the  conduit.  The  conduits  were  assumed  to  be  5  x  5  m,  either  square 
or  semi-elliptical,  and  renting  directly  on  the  incompressible  foundation.  By  exaggera¬ 
ting  conventional  shapes  of  conduits,  a  better  insight  into  the  effect  of  the  geometry 
of  the  conduit  was  obtained. 


2.2  CONDUIT  IN  HOMOGENEOUS  FILL 

Three  cases  with  homogeneous  fill,  illustrated  in  Fig.  1,  were  mvestigated.  The 
fill  is  assumed  to  have  the  following  properties: 

Unit  weight  =  1  900  kg/cu  m 

Young’s  modulus  =  200  kg/sq  cm 

Poisson’s  ratio  =  0.35 

For  the  sake  of  analysis,  vertical,  rigid-slippery  boundaries  are  assumed  20  m  from  the 
centerline  of  the  conduit. 

Case  A  -  Semi- elliptical  Conduit  with  Ratio  of  Height  of  Embankment  to  Width  of 
Conduit  H/D  =  5 

Fig.  2(a)  shows  (1)  the  pertinent  dimensions  and  (2)  the  principal  stresses  in 
the  embankment.  With  the  height  of  the  embankment  H  =  25  m  and  the  width  of  the 
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conduit  D  =  5  m,  the  ratio  is  H/D  =  5.  As  can  be  seen,  no  tensile  stresses  are 
induced  in  ihe  fill  adjacent  to  this  conduit,  but  there  is  a  zone  of  tensile  stresses 
along  the  surface  of  the  fill.  This  zone  extends  about  10  meters  on  both  sides  of 
the  centerline  of  the  conduit.  Its  maximum  depth  is  1.2  maters  and  the  maximum 

lensile  stress  is  only  0.1  kg/sq  cm.  Both  occur  along  the  centerline  of  the  conduit. 
By  comparison  with  the  investigations  by  Covarrubtas  (1969),  it  is  concluded  that 
tensile  stresses  of  this  order  of  magnitude  normally  will  not  produce  tensile  cracks. 

In  Fig.  2(b)  are  plotted  the  normal  stresses  on  the  surface  of  the  conduit, 
which  have  values  of  pv  =  6.6  kg/sq  cm  at  the  center  of  the  crown  and  p^  =  2.7  kg/sq  cm 
at  mid-height  of  the  conduit.  These  values  define  the  loading  ratio  Py/p^  =  2.4. 


Case  B  —  Square  Conduit  with  H/D  =  5 

Fig.  3(a)  shows  (1)  the  pertinent  dimensions  and  (2)  the  principal  stresses  in 
the  embankment.  Adjacent  to  the  upper  portion  of  the  sides  of  the  conduit  there 
are  small  tension  zones  with  a  maximum  tensile  stress  of  about  0.9  kg/sq  cm  along 
the  upper  edges.  The  tension  zone  along  the  surface  of  the  fill  is  similar  to  that 
in  the  preceding  case,  with  a  maxim  m  tensile  stress  at  the  surface  of  about 
0.2  kg/sq  cm  along  the  centerline  of  the  conduit.  Depending  chiefly  on  the  stress- 
strain  characteristics  of  the  fill,  the  tension  zone  adjacent  to  the  side  walls  of  the 
conduit  might  induce  “piping”  along  the  conduit. 

In  Fig.  3(b)  are  plotted  the  normal  stresses  on  the  surface  of  the  conduit. 

On  the  roof  the  normal  stress  ranges  from  about  5  kg/sq  cm  along  the  centerline 
to  about  8  kg/sq  cm  along  the  edges,  and  along  the  middle  of  the  sides  it  is  about 
2  kg/sq  cm.  The  loading  ratio  pv/ph  =  5/2  =  2.5  is  slightly  greater  than  for  the 
elliptical  conduit,  Case  A. 

Case  C  -  Square  Conduit  with  H/D  =  10 

To  investigate  the  effect  of  a  much  larger  H/D  ratio  on  the  tension  zone 
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adjacent  to  the  square  conduit  without  haying  to  resort  to  preparation  of  a  new 
set  of  data  cards  for  the  computer  program,  a  surcharge  of  4.75  kg/sq  cm.  which 
is  equivalent  to  an  additional  25 -meter  height  of  fill,  was  superposed  on  Case  B. 

For  practical  purposes  this  procedure  yields  the  same  stresses  in  the  immediate  vi¬ 
cinity  of  the  conduit  as  would  be  obtained  for  a  50- meter  high  fill  to  start  witn. 

Fig.  4(a)  shows  the  distribution  of  principal  stresses  in  the  vicinity  of  the 
conduit.  It  can  be  seen  that  the  tension  zones  adjacent  to  the  upper  portion  of 
the  sides  are  smaller  than  for  Case  B,  but  the  maximum  tensile  stress  of  1.8  kg/sq  cm 
is  twice  as  great  as  for  Case  B.  Although  the  much  greater  maximum  tensile  stress 
suggests  a  more  dangerous  condition,  the  greater  height  of  fill  would  also  prodr  se 
a  greater  tendency  for  the  crack  to  close  by  creep  deformation  in  the  fill. 

In  Fig.  4(b)  .*re  plotted  the  normal  stresses  on  the  surface  of  the  conduit. 
On  the  roof,  the  normal  stress  ranges  from  about  1 1  kg/sq  cm  along  the  center- 
line  to  about  15  kg/sq  cm  along  the  edges,  and  along  the  middle  of  the  sides  it 
is  about  4.2  kg/sq  cm.  The  loading  ratio  pv/pjj  =  2.5  is  the  same  as  for  the  pre¬ 
ceding  Case  B. 


2.3  USE  OF  A  COMPRESSIBLE  ZONE  TO  CONTROL  LOADING  ON  CONDUIT 

The  most  effective  method  for  improving  the  stress  distfibution  around  a 
conduit,  i.e.  to  achieve  a  loading  ratio  Py/p^  close  to  one,  is  the  addition  of  a 
more  compressible  zone  on  top  of  the  conduit.  Three  cases  have  been  investigated 
in  which  the  dimensions  of  the  compressible  zone  are  assumed  to  be  identical  with 
the  dimensions  of  the  5  x  5- meter  square  conduit,  and  using  compressibility  ratios 
Ec/Ef,  i.e.  ratios  of  the  modulus  of  elasticity  Ec  of  the  more  compressible  zone 
to  the  modulus  of  elasticity  Ef  of  the  fill  material,  of  0.5,  0.2  and  0.1.  The  unit 
weight  and  Poisson’s  ratio  of  the  compressible  zone  are  assumed  equal  to  those  of 
the  fill  material. 
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Case  D  —  Compressibility  Ratio  Eg/Ef  =  0.5 

Fig.  5(a)  shows  (I)  the  pertinent  dimensions  and  (2)  the  principal  stresses 
in  the  fill.  Adjacent  to  the  upper  portions  of  the  sides  of  the  conduit  there  are 
small  tension  zones  with  a  maximum  tensile  stress  of  about  1.2  kg/sq  cm  along 
the  upper  edges.  The  tension  zone  along  the  ciest  of  the  fill  is  similar  to  the 
case  of  the  homogeneous  fill  in  Fig  3(a). 

As  can  be  seen  in  Fig.  5(b)  the  vertical  stresses  on  the  roof  of  the  con¬ 
duit  increase  from  a  value  of  3.8  kg/sq  cm  along  the  centerline  to  about  5.G  kg/sq  cm 
along  the  edges.  The  maximum  horizontal  stress  along  the  sides,  of  about 
2.1  kg/sq  cm  occurs  at  about  mid-height  on  the  conduit.  The  loading  ratio,  pv/pjj  =  1.8, 
is  about  70°/o  of  the  corresponding  value  in  Case  B  without  the  compressible  zone. 


Case  E  -  Compressibility  Raxio  Eg/Ef  -  0.2 

For  the  same  dimensions  as  in  Fig.  5(a),  but  for  a  compressibility  ratio  of 
0.2,  the  distribution  of  the  principal  stresses  and  of  the  tension  zones  along  the 
crest  of  the  fill  and  along  the  sides  of  the  conduit  are  shown  in  Fig.  6(a).  Now 
there  is  no  tension  zone  along  the  middle  of  the  crest,  but  there  are  tension  zones 
along  both  ends  of  the  crest,  with  maximum  tensile  stresses  of  about  U.06  kg/sq  cm 
at  the  ends.  The  tension  zones  along  the  sides  of  the  conduit  have  a  maximum  tensile 
stress  of  about  1.3  kg/sq  cm  along  the  upper  edges. 

Fig.  6(b)  shows  the  distribution  of  normal  stresses  on  the  surface  of  the  con¬ 
duit.  On  the  roof  the  normal  stresses  range  from  2.3  kg/sq  cm  along  the  centerline 
to  about  2.9  kg/sq  cm  along  the  edges.  The  maximum  horizontal  stress  along  the 
sides  is  about  2.3  kg/sq  cm  and  occurs  at  about  mid-height  on  the  conduit.  There¬ 
fore,  one  achieves  in  this  case  the  ideal  loading  ratio  Py/p^  =  1.0. 

Case  F  —  Compressibility  Ratio  Ep/Ef  =  0.1 

For  the  same  dimensions  as  in  Fig.  5(a),  but  for  a  compressibility  latio  of 
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0.1,  the  distr.bution  of  the  principal  stresses  and  of  the  tension  zones  along  the  crest 
of  the  fill  and  along  the  sides  of  the  conduit  are  shown  in  Fig.  7(a),  Again  there 
are  two  tension  zones  along  both  ends  of  the  crest,  with  a  maximum  tensile  stress 
of  about  0.1  kg/sq  cm  at  the  ends.  The  maximum  tensile  stress  along  the  sides  of 
the  conduit  is  about  1.3  kg/sq  cm  and  is  again  located  along  the  upper  edges. 

Fig  7(b)  shows  the  distribution  of  noimal  stresses  on  the  surface  of  the  con¬ 
duit.  On  the  roof  the  normal  stresses  range  from  1.8  kg/sq  cm  along  the  center- 
line  to  2,0  kg/sq  cm  along  the  edges,  and  the  maximum  horizontal  stress  along 
the  sides  of  about  2.5  kg/sq  cm  occurs  about  mid-height  of  the  conduit.  Thus, 
the  loading  ratio  is  pv/p}j  =  0.7,  i.e.  less  than  unity. 


2.4  DISCUSSION  OF  RESULTS 

Table  1  is  a  summary  of  the  most  important  results  in  this  chapter. 

Tension  Zones  Along  Crest  of  Dam  —  For  the  investigated  H/D  ratios  of  5  and  10 
the  maximum  induced  tensile  stresses  axe  so  small  that  they  would  not  produce  ten¬ 
sion  cracks. 

Tension  Zones  Adjacent  to  Conduit  —  For  square  conduits  small  tension  zones  occur 
in  the  fill  adjacent  to  the  upper  portions  of  the  sides  of  the  conduits.  The  tensile 

stresses  that  develop  may  cause  formation  of  tension  cracks  in  the  fill.  Such  cracks, 

combined  with  the  possibility  that  the  soil  will  pull  away  from  the  side  walls  near 
the  upper  edges,  may  induce  “piping"  along  the  conduit  walls. 

With  increasing  ratio  H/D  the  tension  zones  adjacent  to  square  conduits  de¬ 
crease  slightly  in  size,  but  the  maximum  tensile  stresr  increases  substantially. 


For  the  semi-elliptical  eonduit  the  analysis  yielded  no  tension  zones  adjacent 
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to  the  conduit.  This  result  suggests  that  sharp  edges  should  be  avoided  on  the  out¬ 
side  of  conduits. 

Effect  of  Highly  Compressible  Zone  over  a  Conduit  —  The  presence  of  a  highly  com¬ 
pressible  zone  over  the  roof  of  a  conduit  effectively  reduces  the  magnitude  of  the 
vertical  stresses  on  the  roof  and  slightly  increases  the  horizontal  stresses  on  the  sides 
of  the  conduit.  By  appropriate  choice  of  the  height  of  the  compressible  zone  and 
of  the  compressibility  ratio  a  desired  ’oading  ratio  can  be  achieved.  For  the  geome¬ 
try  illustrated  in  Fig.  1,  the  results  plotted  in  Fig.  8  show  tuat  for  a  compressibility 
ratio  Ec/Ef  of  about  0.2,  one  achieves  a  loading  ratio  pvlp^  -  1.0,  i.e.  a  hydrostatic 
loading  condition. 
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CHAPTER  3 

TENSION  ZONES  IN  AN  EMBANKMENT  AND  ITS  FOUNDATION  CAUSED  BY 
A  RIGIDLY  SUPPORTED  CUTOFF  WALL 


3.1  GENERAL 

In  this  chapter  are  presented  the  results  of  investigations  of  the  following  four 

cases: 

Case  A.  Stress  distribution  in  an  embankment  and  in  its  compressible  foun¬ 
dation,  with  a  thick  rigid-rough  cutoff  wall  through  the  foundation,  i.e.  with 
zero  relative  displacements  between  the  wall  and  the  adjacent  soil. 

CaT'e  B.  Same  as  Case  A,  except  that  the  surfaces  of  the  cutoff  wall  are  as¬ 
sumed  to  be  “slippery”,  i.e.  no  friction  forces  are  transmitted  between  the 
wall  and  the  soil. 

Case  C.  Same  as  Case  A,  except  that  a  zone  of  highly  compressible  material 
is  placed  immediately  above  the  top  of  the  cutoff  wall. 

Case  D,  Same  as  Case  A,  except  that  the  rigid-rough  cutoff  wall  is  assumed 
to  have  zero  thickness. 

The  dist  ibution  of  stresses  acting  on  the  sides  of  the  cutoff  wall  is  studied  in 
an  effort  to  determine  the  order  of  magnitude  of  the  vertical  stresses  in  the  wall 
which  are  produced  by  the  weight  of  the  embankment. 

For  Cases  A,  B  and  C,  the  geometry  f  the  embankment  is  shown  in  Fig.  9(a). 
It  is  100  m  high,  with  side  slopes  of  1  (vertical)  on  3  (horizontal),  and  rests  on  a 
100-m-thick  foundation  stratum  which  has  the  same  compressibility  as  the  embank¬ 
ment  and  which  in  turn  is  underlain  by  a  horizontal  rigid-rough  surface.  Along  the 
centerline  a  10-m-wide  rigid  cutoff  wall  extends  through  the  compressible  foundation 
stratum,  i.e.  from  the  base  of  the  dam  down  to  the  rigid  underlying  rock. 

For  Case  D,  the  geometry  of  the  embankment  is  shown  in  Fig.  9(bj.  It  has 
the  same  height,  side  slopes  and  foundation  stratum  as  in  cases  A,  B  and  (',  but  the 
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rigid-rough  wall  in  assumed  to  have  a  zero  thickness. 

(-'or  convenience  of  analysis,  in  all  eases  vertical  “slippery”  boundaries  are 
assumed  on  both  sides  of  tlie  embankment  at  distances  of  500  m  from  the  center- 
line.  In  addition,  it  is  assumed  that  the  compressible  foundation  deforms  only  under 
the  weight  of  the  embankment  and  not  under  its  own  weight,  which  is  equivalent 
to  assume  zero  unit  weight  for  the  foundation  soil.  The  assumed  material  proper¬ 
ties  are  listed  below: 


Zone 

Young’s 

Modulus 

Poisson’s 

Unit 

weight 

kg/sq  cm 

ratio 

kg/cu  m 

Embankment 

200 

0.35 

1  900 

Foundation 

200 

0.35 

0 

32  CASE  A.  THICK  CUTOFF  WALL  WITH  ROUGH  SIDES 

In  Case  A  the  wall  is  assumed  to  be  “rough,”  i.e.  no  relative  movements 
can  develop  between  the  wall  and  the  adjacent  soil. 

Fig.  10(a)  shows  the  displacements  of  points  within  the  embankment  and 
within  the  foundation.  The  maximum  settlement  of  5.72  m  occurs  at  the  center- 
line  of  the  crest.  The  maximum  settlement  of  the  foundation  surface  (base  of 
embankment)  of  3.83  m  occurs  about  90  m  from  the  centerline. 

Fig.  10(b)  shows  t’-e  distribution  of  the  principal  stresses,  and  the  locations 
of  the  tension  zones  in  the  embankment  and  in  the  foundation.  The  tension  zone 
m  the  top  of  the  embankment  is  roughly  triangular  and  extends  from  the  crest  to 
about  72  m  down  the  slopes  and  to  a  depth  of  about  50  m  beneath  the  crest. 
The  maximum  tensile  stress  of  about  2.4  kg/sq  cm  occurs  on  the  slope  at  a  dis¬ 
tance  of  approximately  20  m  from  the  crest. 

The  tension  zones  in  the  foundation  stratum  occur  adjacent  to  the  upper  part 
of  the  cutoff  wall.  Thc\  extend  over  a  depth  of  40  m,  with  a  maximum  width 
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of  5  3  m.  The  maximum  tensile  stress  of  about  13  kg/sq  cm  occurs  adjacent  to 
the  top  edges  of  the  cutoff  wall. 

The  vertr  .tresses  oy  on  the  top  surface  of  the  wall  are  about  44  kg/sq  cm, 
which  is  about  2.3  times  larger  than  the  vertical  stress  under  a  100-m-high  column 
of  embankment  material. 

In  Fig.  10(c)  are  shown  (1)  the  distribution  of  effective  horizontal  stresses 
acting  on  the  sides  of  the  wall,  (2)  the  potential  maximum  shearing  resistance  aiong 
the  wall,  and  (3)  the  computed  shear  stresses  along  the  wall  produced  by  the  weight 
of  the  embankment.  The  stresses  (1)  are  the  sum  of:  (a)  the  in  situ  effective  hor¬ 
izontal  stresses  in  the  foundation  stratum  (computed  for  a  submerged  specific  gravity 
of  unity  and  a  coefficient  of  earth  pressure  at  rest  K0  =  0.5);  and  (b)  the  horizon¬ 
tal  stresses  p.oduced  by  the  weight  of  the  embankment.  The  stresses  (2)  are  com¬ 
puted  assuming  a  coefficient  of  wall  friction  of  0.5.  The  shear  stresses  (3)  are  com¬ 
puted  by  means  of  the  finite  element  method. 

The  effective  horizontal  stresses  acting  on  the  sides  of  the  wall  are  compressive 
over  almost  the  entire  length  of  the  wall,  but  near  the  top  of  the  wall,  they  decrease 
abruptly  and  rapidly  and  become  tensile  stresses  over  the  upper  5  m  of  the  wall.  The 
potential  maximum  shearing  resistance  is  sensibly  uniform  over  most  of  the  height  of 
the  cutoff  wall,  with  a  maximum  value  of  about  4.3  kg/sq  cm  approximately  at  mid¬ 
height.  The  shear  stresses  produced  by  the  weight  of  the  embankment  exceed  21  kg/sq  cm 
at  the  top  of  the  wall,  and  they  decrease  rapidly  with  depth.  At  a  depth  of  45  m 
below  the  base  of  the  embankment  both  the  shear  stress  and  the  potential  max  num 
shearing  resistance  are  equal,  and  below  that  depth  the  shear  stresses  are  smaller  than 
the  potential  maximum  shearing  resistance. 

Assuming  that  the  shear  stresses  which  are  transfered  from  the  surrounding  soil 
mass  to  the  wall  are  equal  to  the  shear  stresses  produced  on  the  wall  by  the  weight 
of  the  embankment,  but  not  greater  than  the  potential  maximum  shearing  resistance, 
and  integrating  those  stresses  over  the  length  of  the  wall,  the  total  force  transmitted 
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to  bedrock  through  the  base  of  the  wall  due  to  negative  skin  friction  over  one  side 
of  the  wall  is  23.700  kg/cm  of  wall  length.  Assuming  a  28-day  compressive  strength 
of  the  concrete  of  200  kg/sq  cm,  one  would  need  a  2.4 -m- thick  wall  to  avoid  crush¬ 
ing  of  the  lowest  portion  of  the  wall  by  the  load  produced  by  negative  skin  friction 
which  would  develop  in  the  course  of  years. 


3.3  CASE  B.  THICK  CUTOFF  WALL  WITH  SLIPPERY  SIDES 

Case  B  is  similar  to  Case  A,  except  that  the  sides  of  the  cutoff  wall  are  as¬ 
sumed  to  be  slippery.  Therefore,  the  shear  stresses  along  the  wall  are  zero  and  the 
soil  adjacent  to  the  wall  can  displace  vertically. 

Fig.  11(a)  shows  the  displacements  of  points  within  the  embankment  and 
within  the  foundation.  The  maximum  settlement  of  6.03  m  occurs  at  the  centerline 
of  the  crest.  The  maximum  settlement  of  the  foundation  surface  (base  of  embank¬ 
ment)  of  3.99  m  oeeurs  about  90  m  from  the  centerline.  Both  of  the  movements 
are  slightly  greater  than  for  Case  A  because  of  the  freedom  of  movement  at  the 
wall. 


Fig.  11(b)  shows  the  distribution  of  the  principal  stresses  and  the  locations 
of  the  tension  zones  in  the  embankment  and  in  the  foundation.  The  tension  zone 
in  the  top  ot  the  embankment  is  roughly  triangular  and  extends  from  the  crest  to 
about  70  m  down  the  slopes  and  to  a  depth  of  bout  52  in  beneath  the  erest.  The 
maximum  tensile  stress  of  about  2  kg/sq  cm  occurs  on  the  slopes  at  a  distance  of 

about  20  m  from  the  crest. 

The  tension  zones  in  the  foundation  stratum  appear  adjacent  to  the  upper 
part  of  the  cutoff  wall,  extending  over  a  depth  of  19  m,  with  a  maximum  width 
of  15  in.  The  maximum  tensile  stress  of  about  30  kg/sq  cm  occurs  adjacent  to 
the  top  edges  of  the  cutoff  wall. 


The  vertical  stresses  on  the  top  surface  of  the  wall  are  about  50  kg/sq  em. 
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i.e.  approximately  2.6  times  larger  than  the  stress  under  a  100~m-high  column  of 
embankment  material.  In  this  case,  since  the  wall  has  been  assumed  to  be  slippery, 
the  wall  is  not  loaded  by  negative  skin  friction. 


3.4  CASE  C.  COMPRESSIBLE  ZONE  ON  TOP  OF  A  THICK  CUTOFF  WALL 
WITH  ROUGH  SIDES 

Case  C  is  similar  to  Case  A,  except  for  the  assumption  of  a  compressible 
zone  on  top  of  the  cutoff  wall  with  Young’s  modulus  10  times  smaller  than  the 
modulus  in  the  embankment  proper,  and  a  10  x  10-m  cross  section. 

Fig.  12(a)  shows  the  displacements  of  points  within  the  embankment  and 
within  the  foundation.  The  maximum  settlement  of  6.13  m  occurs  at  the  centerline 
of  the  crest.  The  maximum  settlement  of  the  foundation  surface  (base  of  embank¬ 
ment)  of  3.90  m  occurs  about  90  m  from  the  centerline. 

Fig.  12(b)  shows  the  distribution  of  principal  stresses  and  the  locations  of 
the  tension  zones  in  the  embankment  and  in  the  foundation.  The  tension  zone  in 
the  top  of  the  embankment  is  almost  triangular  and  extends  from  the  crest  about 
70  m  along  the  slopes  and  to  a  depth  of  about  45  m  beneath  the  crest.  The 
maximum  tensile  stress  of  about  1.5  kg/sq  cm  occurs  on  the  slopes  at  a  distance 
of  approximately  20  m  from  the  crest. 

The  tension  zones  in  the  foundation  occur  adjacent  to  the  upper  half  of  the 

cutoff  wall  and  have  a  maximum  width  of  13  m.  The  maximum  tensile  stress  of 

about  16  kg/sq  cm  occurs  adjacent  to  the  upper  edges  of  the  cutoff  wall. 

The  vertical  stresses  on  the  top  surface  of  the  wall  are  about  8.5  kg/sq  cm, 

which  is  only  about  one-half  of  the  vertical  stress  under  a  100-m-high  column  of 
embankment  material. 
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In  f  ig.  12(c)  arc  shown’  (1)  the  distribution  of  effective  stresses  acting  on 

the  sides  of  the  wall,  with  the  same  assumptions  of  specific  gravity  and  KQ  as  in 

Case  A,  (2)  the  potential  maximum  shearing  resistance  along  the  wall  using  the  same 
value  of  the  coefficient  of  wall  friction  as  in  Case  A,  and  (3)  the  shear  stresses  along 
the  wall  caused  b>  the  weight  of  the  embankment. 

The  effective  horizontal  stresses,  which  are  compressible  over  almost  the  entire 
height  of  the  wall,  decrease  abruptly  near  the  top  and  become  tensile  stresses  over  the 
top  7  m  of  the  wall.  The  potential  maximum  shearing  resistance  has  two  maxima: 

(1)  about  5.3  kg/sq  cm  approximately  15  m  below  the  top  of  the  wall,  and  (2)  about 
4.5  kg/sq  cm  approximately  70  m  beiow  the  top  of  the  wall.  At  the  top  of  the  wall, 

the  shear  stresses  caused  by  the  weight  of  the  embankment  exceed  28  kg/sq  cm  and 

decrease  rapidly  with  depth.  At  a  depth  of  48  m  below  the  top  of  the  wall,  both 
the  shear  stress  ar.d  the  shearing  resistance  are  equal.  Below  that  depth  the  shear  stresses 
are  smaller  than  the  potential  maximum  shearing  resistance.  The  total  force  transmitted 
to  the  base  of  the  wall  by  negative  skin  friction  is  26,520  kg/cm  for  each  side  of  the 
wall. 


As  compared  to  Case  A,  for  Case  C  (1)  the  computed  settlements  are  larger, 

(2)  the  tension  zone  in  the  embankment  is  smaller,  (3)  the  maximum  tensile  stress  in 
the  embankment  is  substantially  smaller,  (4)  the  tension  zone  in  the  foundation  stratum 
is  slightly  deeper,  and  (5)  the  maximum  tensile  stress  in  the  foundation  is  greater. 


3.5  CASE  D.  THIN  CUTOFF  WALL  WITH  ROUGH  SIDES 

Case  D  is  similar  to  Case  A,  except  that  the  cutoff  wall  which  extends  through 
the  foundation  is  assumed  to  have  zero  thickness. 

Fig.  13(a)  shows  the  displacements  of  points  within  the  embankment  and  within 
the  foundation.  The  maximum  settlement  of  6.37  m  occurs  at  the  centerline  of  the 
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crest.  The  maximum  settlement  of  the  foundation  surface  (base  of  embankment) 
of  3.98  m  occurs  about  90  m  from  the  centerline. 

Fig.  13(b)  shows  the  distribution  of  principal  stresses  and  the  locations  of  the 
tension  zones  in  the  embankment  and  in  the  foundation.  The  tension  zone  in  the 
top  of  the  embankment  is  roughly  triangular  and  extends  from  the  crest  about  60  m 
along  the  slopes  and  to  a  depth  of  about  50  m  beneath  the  crest.  The  maximum  ten¬ 
sile  stress  of  about  1.5  kg/sq  cm  occurs  on  the  slopes  approximately  20  m  from  the 
crest. 


The  tension  zones  in  the  foundation  extend  adjacent  to  the  upper  44  m  of  the 
cutoff  wall  and  have  a  maximum  width  of  13  m.  The  maximum  tensile  stress  of 
about  15  kg/sq  cm  occurs  adjacent  to  the  upper  edges  of  the  cutoff  wall. 

In  Fig.  13(c)  are  shewn:  (1)  the  distribution  of  effective  stresses  acting  on  the 
sides  of  the  wall,  with  the  same  assumptions  of  specific  gravity  and  K0  as  in  Case  A, 
(2)  the  potential  maximum  shearing  resistance  along  the  wall  with  the  same  value  of 
the  coefficient  of  wall  friction  as  in  Case  A,  and  (3)  the  shear  stresses  along  the  wall 
caused  by  the  weight  of  the  embankment. 

The  effective  horizontal  stresses  acting  on  the  sides  of  the  wall,  which  are 
compressive  over  almost  the  entire  length  of  the  wall,  decrease  abruptly  near  the 
top  and  become  tensile  stresses  over  the  top  3  m  of  the  wall.  The  potential  max¬ 
imum  shearing  resistance  is  sensibly  uniform  over  almost  the  entire  height  of  the 
cutoff  wall,  with  a  maximum  value  of  about  4.5  kg/sq  cm  approximately  65  m  be¬ 
low  the  top  of  the  wall.  The  shear  stresses  caused  by  the  weight  of  the  embank¬ 
ment  exceed  21  kg/sq  cm  at  the  top  of  the  wall,  and  they  decrease  rapidly  with  depth. 
At  a  depth  of  45  m  below  the  top  of  the  wall,  both  the  shear  stress  and  the  shearing 
resistance  are  equal,  and  below  that  depth  the  shear  stresses  are  smaller  than  the 
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potential  maximum  shearing  resistance.  With  the  same  assumptions  as  stated  in  Case 
A,  the  total  force  transmitted  to  the  foundation  through  the  base  of  the  wall,  caused 
hy  negative  skin  friction,  is  25,400  kg/cm  for  each  side  of  the  wall. 

3.6  DISCUSSION  OF  RESULTS 

Table  2  is  a  summary  of  the  important  numerical  results  in  this  chapter.  For 
the  geometry  studied,  which  is  shown  in  Figs.  9(a)  and  9(b),  the  most  important  con¬ 
clusions  are: 

(1)  The  presence  of  a  rigid  cutoff  wall  induces  substantial  tension  zones  in 
the  embankment  and  in  its  foundation. 

(2)  Large  forces  can  be  transmitted  by  negative  skin  friction  from  the  soil 
to  the  cutoff  wall  which  may  exceed  the  strength  of  the  wall  in  its  lower 
portion. 

Tension  Zones  in  the  Embankment  —  The  tension  zone  produced  in  the  em¬ 
bankment  by  a  rigid  cutoff  wall  is  always  located  at  the  top  or  the  embankment, 
with  the  maximum  tensile  stresses  along  the  slopes. 

The  largest  tensile  stress  in  the  embankment  occurs  for  the  case  of  a  thick 
cutoff  wall  with  rigid-rough  sides,  Case  A,  while  the  smallest  tensile  stresses  occur 
for  Case  C,  a  thick  cutoff  wall  with  rigid-rough  sides  and  with  a  highly  compressible 
zone  on  top  of  the  wall,  and  for  Case  D,  a  thin  cutoff  wall  with  rigid-rough  sides. 
The  computed  values  of  the  maximum  tensile  stress  ranged  from  1.5  kg/sq  cm  (Cases 
C  and  D)  to  2.4  kg/sq  cm  (Case  A),  which  are  relatively  small  as  compared  with 
those  computed  by  Covarrubias  (1969)  in  actual  dams  that  have  cracked. 


Tension  Zones  in  the  Foundation  Stratum  —  The  tension  zones  in  the 
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foundation  stratum  are  located  adjacent  to  the  upper  part  of  the  cutoff  wall,  with 
the  maximum  tensile  stress  occurring  along  the  top  edges  of  the  wall.  The  magnitude 
of  the  maximum  tensile  stress  is  greatest  for  the  case  of  a  thick  cutoff  wall  with  rigid- 
slippery  sides,  Case  B.  There  are  no  large  differences  among  the  other  cases.  The 
maximum  tensile  stresses  range  from  13  to  30  kg/sq  cm,  values  which  no  soil  can  with¬ 
stand.  However,  at  such  depths  the  plasticity  and  time-dependent  properties  of  soils 
may  play  an  important  role  in  effecting  substantial  redistribution  of  stresses. 

Vertical  Stresses  on  the  Top  of  the  Wall  —  The  vertical  stresses  acting  on  the 
top  surface  of  the  wall  are  effectively  reduced  by  the  presence  of  a  zone  of  more  com¬ 
pressible  material  over  the  top  of  the  wall.  When  the  more  compressible  zone  is  ten 
times  more  compiessible  than  the  surrounding  fill,  Case  C,  the  vertical  stress  is  only  25°/D 
of  the  corresponding  stress  in  a  homogeneous  embankment,  Case  A,  or  0.45  times  the 
weight  of  a  column  of  soil  equal  to  the  height  of  the  embankment. 

Stresses  and  Friction  on  the  Side  of  the  Wall  —  The  effective  horizontal  stresses 
that  act  against  the  wall  are  tensile  stresses  over  a  small  region  along  the  upper  part  of 
the  wall  and  they  change  rapidly  with  depth  to  become  compressive  at  3  to  8-m  depth, 
in  homogeneous  embankments  (Cases  A,  B  and  D)  the  distribution  of  the  effective  hor¬ 
izontal  stresses  is  rather  uniform  over  most  of  the  wall,  while  for  Case  C  it  shows  two 
maxima  approximately  at  the  third  points  of  the  wall. 

In  all  cases,  the  computed  shear  stresses  on  the  side  of  the  wall  caused  by 
the  weight  of  the  embankment,  are  extremely  large  along  the  upper  part  of  the  cutoff 
wall  and  decrease  rapidly  with  depth. 

The  negative  skin  friction  is  slightly  larger  for  the  case  in  which  a  highly  com¬ 
pressible  zone  is  assumed  on  top  of  the  wall,  Case  C.  However,  the  presence  of  the 
compressible  zone  results  in  a  small  downward  force  on  top  of  the  wall.  The  dif¬ 
ference  in  negative  skin  friction  between  Case  C  and  the  other  cases  is  not  sub¬ 
stantial. 


28 


TABLE  1 

RESULTS  OF  FINITE  ELEMENT  ANALYSIS  OF  EMBANKMENT 
SURROUNDING  RIGiD  CONDUITS 


Ca si 

H/D 

Compress¬ 

ibility 

ratio 

Ec/Ef 

Maximum  tensile 

stress 

in  the  crest 

Meximum  tensile 
stress  adjacent 
to  conduit 

Loading 

ratio 

Pv/Ph 

kg/sq  cm 

kg/sq  cm 

A 

5 

1 

0.1 

- 

2.4 

B 

5 

1 

0.2 

0.9 

2.5 

C 

10 

1 

- 

1.8 

2.5 

D 

5 

0.5 

0.05 

1.2 

1.8 

E 

5 

0.2 

0.06 

1.3 

1.0 

F 

5 

0.1 

0.1 

1.3 

0.7 

TABLE  2 

RESULTS  OF  FINITE  ELEMENT  ANALYSIS  OF  EMBANKMENT  WITH 
RIGIDLY  SUPPORTED  CUTOFF  WALL 


Case 

Tension  zone  in 
embankment 

Tension  zon.  in 
foundation 

1 

sr.iax 

Tmax 

F 

Maximum 

settlement 

width 

depth 

maximum 

tensile 

stress 

width 

depth 

maximum 

tensile 

stress 

west 

foun¬ 

dation 

surface 

m 

m 

kg/sq  cm 

m 

m 

kg/sq  cm 

kg/sq  cm 

kg/cm 

m 

m 

A 

72 

50 

2.4 

13 

40 

13 

44 

4.3 

21 

23760 

5.72 

3.83 

B 

70 

52 

2 

15 

19 

30 

50 

- 

- 

- 

6.03 

3.97 

C 

70 

45 

1.5 

13 

50 

16 

8.5 

5.3 

28 

26520 

6.13 

3.90 

D 

60 

50 

1.5 

13 

44 

15 

- 

4.5 

l_£L_ 

25400 

6.37 

3.98 

CASE  A- Semi-elliptical  conduit,  H/D  =  5 
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Fig  2(a)  Principal  stresses  in  fill  surrounding  semi-elliptical  conduit 


Fig.  2(b)  Normal  stresses  acting  on  the  sides  of  semi-elliptical  conduit 
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Dimensions  ore  in  meters 
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Fig.5(a)  Principal  stresses  in  fill  surrounding  square  conduit;  Compressibility  ratio  of  1/2 
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Fig.  5(b)  Normal  stresses  acting  on  the  sides  of  square  conduit;  Compressibility  ratio  of  1/2 
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Fig.  6(a)  Principal  stresses  in  fill  surrounding  square  conduit;  Compressibility  ratio  of  1/5 
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big.  6(b)  Normal  stresses  acting  on  the  sides  of  square  conduit;  Compressibility  ratio  of  1/5 
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Fig.  7(a)  Principal  stresses  in  fill  surrounding  square  conduit;  Compressibility  ratio  of  1/10 


Fig.  7(  b)  Normal  stresses  acting  on  the  sides  of  square  conduit;  Compressibility  ratio  of  1/10 
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fig.  &  Loading  ratio  for  square  conduit  as  function  of  compressibility  ratio 
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